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ABSTRACT 

We report the discovery of TeV gamma-ray emission coincident with the shell-type radio supernova 
remnant (SNR) CTA 1 using the VERITAS gamma-ray observatory. The source, VER J0006+729, 
was detected as a 6.5 standard deviation excess over background and shows an extended morphology, 
approximated by a two-dimensional Gaussian of semi-major (semi-minor) axis 0°.30 (0°.24) and a 
centroid 5' from the Fermi gamma-ray pulsar PSR J0007-I-7303 and its X-ray pulsar wind nebula 
(PWN). The photon spectrum is well described by a power-law dN/dE = Nq[E/2> TeV)^'", with a 
differential spectral index of F = 2.2±0.2stat±0.3sys, and normalization Nq = (9.1 ± 1.3stat ± l-7sys) x 
10"^'' cm~^ s~^ TeV~^. The integral flux, — 4.0 x 10~^^ erg cm~^ s~^ above 1 TeV, corresponds 
to 0.2% of the pulsar spin-down power at 1.4 kpc. The energetics, co-location with the SNR, and the 
relatively small extent of the TeV emission strongly argue for the PWN origin of the TeV photons. 
We consider the origin of the TeV emission in CTA 1. 

Subject headings: supernova remnants: individual (SNR G119.5+10.2 = VER J0006+729) — gamma- 
rays: individual (3EG JOOlO+7309), pulsars: individual (PSR J0007-K7303), indi- 
vidual (RX J0007.0+7303) 
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1. INTRODUCTION 

There are many possible associations of gamma-ray 
sources with supernova remnants (SNRs) . These gamma 
rays could come from shock acceleration in the shell, 
a pulsar associated with the SNR, or a pulsar wind 
nebula (PWN) surrounding the pulsar. Such gamma- 
ray/SNR associations date back to COS B observa- 
tions of SNRs coin cident with OB stellar associations 
(iMontmerle I Il979f ). Observations of the Galaxy by 
EGRET in the energy range 30 MeV-30 GeV revealed 
~ 19 unidentified sources at low Galactic latitudes that 
were found t o be spatially correlated with mostly shell- 
type SNRs CTorre s et alTl200l . One such source was 
3EG JOOlO +7309, with a relativ ely smah 95% error cir- 
cle of 28' (|Hartman et alJ [19991 ) . that was spatially co- 
incident with the SNR CTA 1 (G119.5-hl0.2) and the 
X-ray point source R X J0007. 0-1-7303, wh ich was postu- 
lated to be a pulsar (i Brazier et al.|[l998l ). The associ- 
ation between 3EG J0010-h7309 and RX J0007.0-K7303 
was found to be plausible, given the lack of flux vari- 
ability seen in 3EG JOOlO-t-7309, its hard spectral index 
(r = 1.58 ± 0.18 between 70 MeV and 2 GeV), and its 
similarity with othe r known pulsars detected by EGRET 
(jBrazier et al.lll998D . 

CTA 1 is a comp osite SNR, discovered by 
[Harris fc Robertsi (jl960( ). with a shell- type struc- 
ture in the radio band and a center-filled morphology at 
X-ray energies. The radio shell is incomplete towards 
the north-west (NW) of the remnant, possibly due to 
rapid expansion of th e shock into a lower-density region 
(|Pineault et al.lll993D . The distance to SNR CTA 1 as 
derived from the asso ciated HI shell is d = 1.4 ± 0.3 kpc 
(jPineault et al.l [19971), the SNR a ge is estimated to be 
~ 1.3 X 10^ yr (iSlane et al.ll2004D. and the diameter o f 
its radio sheh is ~ 1°.8 (jSieber. Salter, fc Maver |[l981[ ). 

Archival X-ray observations of SNR CTA 1 in the 
5-10 keV band show non-thermal diffuse emission of 
low surface brightness in the center of the remnant, 
likely correspondin g to a PWN dr iven by an ac- 
tive neutron star ()Slane et al.l I1997D . The neutron 
star candidate RX J0007. 0-1-7303 is a faint source lo- 
cated at th e brightest part of the synchrotron emission 
(jSeward. S chmidt & Slanc 1995). A Ghandra image of 
RX J0007. 0-1-7303 provides further evidence that this 
source is an energetic rotation-powered pulsar, resolving 
a central point sourc e, a compact nebula, and a bent jet 
(jHalpern et al.ll2004[ ). An initial observation with XMM- 
Newton in 2002 found the X-ray spectrum of the central 
source to be consistent with that of a neutron star, a l- 
though no pulsations were detected (Sl ane et al.ll2004D . 
Based on these initial X-ray observations, the spin-down 
luminosity of the underlying pulsar was estimated to be 
in the range 10'^^ — 10'^^ ergs s~^, supporting the iden- 
tificati on of the EG RET source 3EG JOOlO+7309 as a 
pulsar (jSlane et al. f 2004; Halpern et al. 2004). 

Eventually, a search for pulsed GeV emission from 
CTA 1 using the data from the Fermi Gamma-Ray Space 
Telescope revealed a highly significant 316 ms signal, 
confirming the origin of 3EG JOOlO-l-7309 ()Abdo et al.l 



Street, Cambridge, MA 02138, USA 

Argonne National Laboratory, 9700 S. Cass Avenue, Argonne, 
IL 60439, USA 



[2001 . The spin-down power was determined to be 
~ 4.5 X 10'^^ erg s""' ^^, which is sufficie nt to power the 
pulsar wind nebula (|Slane et al.l 120041 ) . Following the 
Fermi discovery of the gamma-ray pulsar, a deep 130 ks 
observation of RX J0007. 0+7303 was carried out with 
XMM-Newton to characterize the timing behavior. The 
X-ray signal of PSR J0007-I-7303 was discovered at a 
statistical significance of A. 7 a in the 0. 5-2 keV band, 
out o f phase with the gam ma-ray pulse (iCaraveo et al.l 
W^- Similar to Geminga (Halpern & Hoh "1992') and 
PSR J1836+592 5 (. Hal pern. Camilo, & Gotthelf 2003 
lAbdo et al.ll20T0l ). PSR J0007-F7303 is also radio quiet. 
GeV emission in the off-pulse p hase interval has also re- 
cently been detected by Fermi (jAbdo et al.ll2012| ). 

Many galactic TeV sources appear to be associated 
with pulsars via their wind nebulae, comprised of rel- 
ativistic wind particles confined by the pressur e of the 
surrounding medium (jCaensler fc Slanel l2006t ). The 
initially highly supersonic wind terminates in a shock 
which can be associated with axisymmetric, toroidal 
struct ures often seen in the X- ray images of PWNe 
(e.g., iKargaltsev fc Pavlovl l2008l ). PWNe now repre- 
sent the most populous c lass of Galactic TeV emitters 
()Hinton fc HofmannllMot) . The non-thermal emission 
seen in PWNe from the radio up to gamma rays be- 
low a few GeV or less is generally interpreted as syn- 
chrotron radiation from the accelerated leptons. The 
emission observed at higher energies, up to several TeV, 
can be produced via inverse Compton (IC) scattering of 
these same high-energy electrons with background pho- 
tons (e.g. the cosmic microwave background (CMB), in- 
frared rad iation from dust, starlight, and synchrotron 
photons) (jAtovan fc Aharonian 1 119961 ). Alternatively, 
hadronic mechanisms may also be responsible for the 
TeV emission, in which case the wind should be com- 
posed of relativistic hadrons that collide with the ambi- 
ent medium and produce pions, with the TeV emission 
coming from tt" decay. To date, however, there has been 
no solid evidence requiring a large contribution to the 
gamma-ray emission from such hadronic processes. 

PWNe experience se veral stages of evolution (e.g. 
iGaensler fc Slanel I2006D . At an early stage the pulsar 
wind freely expands into the SN ejecta. For a slowly 
moving pulsar the PWN is approximately centered on 
the pulsar while for a supersonically moving pulsar the 
PWN will take a cometary shape. At later times, the 
PWN is compressed by the reverse SNR shock and may 
be displaced significantly from the pulsar if the reverse 
shock is asymmetric. Such crushed and displaced PWNe 
have been dubbed relic PWNe. 

The X-ray and gamma-ray observations of CTA 1 
suggest that the extended non-thermal emission around 
the gamma-ray pulsar is a synchrotron PWN. Mo- 
tivated by these ob s ervati ons, model calculations by 
iZhang. Jiang fc Lin I (j2009t ) suggested that the TeV 
emission is largely produced by the PWN, and that the 
level of emission should be detectable at TeV energies 
by VERITAS. For IC scatt ering off the PWN electrons, 
IZhang. Jiang fc Lin I H2OO90 predicted a gamma-ray flux 
F^il - 30 TeV) ~ 1.1 X lO'^^ gj.g cmr'^ s^^. Previ- 
ous TeV observations of CTA 1 by the earlier imaging 
atmospheric Cherenkov telescopes gave upper limits, as 



follows: 2.64 x 10 photons cm 



[E > 250 GeV) 
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by CAT ()Khelifi et al.|[200l . 1.25 x 10'" photons cm-^ 
s-i {E > 620 GeV) by Whipple (|Hall et aLII200l . and 
1.09 X 1 0"^^ photons cm- ^ s'^ {E > 1.3 TeV) by 
HEGRA (IRowell et al.ll200l . 

In this paper, we report the VERITAS detection of 
TeV emission from the central region of CTA 1. 

2. VERITAS INSTRUMENT & OBSERVATIONS 

The Very Energetic Radiation Imaging Telescope Ar- 
ray System (VERITAS) uses ground-based detection 
techniques pio neered by its prede cessor, the Whipple 
10m Telescope (jWeekes et al.1119891) . to explore the Uni- 
verse in very high-energy (VHE) gamma rays from ^ 
100 GeV to ~ 30 TeV. The VERITAS telescope ar- 
ray consists of four 12-m diameter Davies-Cotton tele- 
scopes and is located at the basecamp of the Fred 
Lawrenc e Whipple Observ atory (FLWO) in southern 
Arizona ()Holder et al.ll20lil ). Flashes of Cherenkov light 
from gamma-ray and cosmic-ray showers are focused by 
a set of mirrors onto a camera located in the focal plane 
of each telescope. Each camera comprises 499 photo- 
multiplier tube pixels and light concentrators arranged 
in a hexagonal pattern with a total field-of-view of 3°. 5. 
Stereoscopic imaging of showers from multiple viewing 
angles allows the determination of the shower core lo- 
cation relative to the array using simple geometric re- 
construction techniques which rely on the fact that the 
major axes of the shower images are projections of the 
shower axis. The combined instrument has an angular 
resolution of < 0°.l (68% containment) and energy reso- 
lution of 15 — 20% for energies > 200 GeV. It can detect 
a source with a flux of 1% of the steady Crab Nebula 
VHE flux at a 5 standard deviation significance level in 
less than 30 hours (Ong et al, 2009). 

A three-level trigger system is used to help eliminate 
background noise. The first trigger occurs at the pixel 
level, requiring the signal to reach a 50 mV threshold 
(corresponding to 4-5 photoelectrons) set by a constant 
fraction discriminator (CFD). The second, telescope- 
level trigger requires at least three adjacent pixels passing 
the CFD trigger to form an image. A third, array-level 
trigger requires simultaneous Cherenkov images in at 
least two telescopes, within a 50 ns time window, which 
then causes a readout of the 500 MSample/s Flash- ADC 
data acquisition system for each pixel. 

CTA 1 was observed over two epochs. The first set 
of observations spanned from September 2010 to Jan- 
uary 2011, with a total livetime of 25 hours 39 min af- 
ter data-quality selection based on weather conditions 
and hardware status. An additional 15 hours 36 min 
of quality-selected data were taken from September to 
Decem ber 2011. Observ ations were taken in "wobble" 
mode (jFomin et al.lll99"3 ). in which the telescope point- 
ing is offset from the source position by some angular 
distance. This method allows for simultaneous collection 
of data and estimation of the background. Due to the ex- 
tended nature of the remnant and expected extension of 
the PWN, an offset distance of 0°.7 was used, larger than 
the typical VERITAS distance of 0°.5. To decrease bias, 
the offset direction was varied between each 20 minute 
run while maintaining the same offset distance, alternat- 
ing between the north, south, east and west directions 
(in the equatorial coordinate system) . Observations were 



taken in a narrow range of zenith angles, 40 — 47°, with 
an average of 42°. 5 for the full dataset. All of the data 
presented here were taken with all four telescopes in the 
array. 

3. ANALYSIS 

The CTA 1 data were processed using st andard VER- 
ITAS analysis techniques, as described in lAcciari et al.l 
(12001 . The cosmic-ray background was suppressed 
efficiently by parametrizing the recorded shower im- 
ages by their principal moments (Hillas 1985), and 
the shower direction and impact parameter were recon- 
structed fror n these images, using stereoscopic meth- 
ods ( see e.g. lAharonian et al.l Il997t iKrawczvnski et ahl 
120061) . Gamma-ray/hadronic shower separation is 
achieved through selection criteria (cuts). Based upo n 
the predicted spectrum of iZhang. Jiang fc Lin I ^OW), 
two sets of standard cuts were used. These cuts were op- 
timized on simulations for sources of ^ 5% of the Crab 
Nebula flux and with moderate and hard spectral in- 
dices (^ -2.5 and -2.0, respectively). For these cuts, at 
least three of the telescopes in the array had to have im- 
ages recorded in the camera, with at least 1200 digital 
counts (^ 240 photo-electrons) for the hard-cut analy- 
sis and 500 digital counts (~ 95 photo-electrons) for the 
moderate-cut analysis. Cuts were also applied to the 
mean scaled length (MSL), mean scaled width (MSW), 
and integrated charge in the signal (size). Finally a 
cut was applied on 9, the angular distance in the field 
of view from the reconstructed arrival direction of the 
event to the putative source location. A cut oi 9 < 0°.09 
(9 < 0°.23) was used for a point-source (extended-source) 
search, with the size of the extended-source cut selected 
a priori. For the analysis presented here, the cuts for the 
moderate- and hard-spectra analysis are MSW < 0.35 
and MSL < 0.7. 

The background was e stimated using the ring back - 
ground model (e.g., see iBerge. Funk fc HintonI [200l . 
with a ring of mean radius 0°.7 and a background to 
source area ratio of 8.0. Regions in the field of view 
containing stars of B magnitude brighter than 6.0 were 
excluded from the background estimation in order to re- 
duce systematic errors. The statistical signi ficance of 
the excess was calculated using Equation 17 of lLi fc Mai 
(1983). The energy threshold for this analysis after ap- 
plying the moderate (hard) cuts is ~ 550 GeV (1 TeV) 
at a zenith angle of 45°, with a systematic error of about 
20% on the energy estimatio n. Two i ndepe nden t analy- 
sis packages, as described bv iCoganI ()2008f) and iDanieTI 
(2008) , were used to reproduce the results presented here 
on CTA 1. 

4. RESULTS 

Figure [T] shows the TeV excess map of the region of 
the sky around CTA 1. The hard-spectrum, extended- 
source analysis produced an excess with a significance of 
7.5 standard deviation (cr) pre-trials, in a search region of 
radius 0°.4 around the pulsar PSR J0007+7303, within 
the radio shell of the SNR CTA 1. Accounting for the two 
sets of applied cuts with two different integration radii, 
and determining the a priori t rials factor by ti ling the 
search region with 0°.04 bins ([Aharonian et al.l [20 06a). 
we estimate the post-trials significance of detection to be 
6.5(7. Overlaid on the TeV image are the high-resolution 
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Figure 1. VERITAS excess map of the region around SNR CTA 1 
using a hard-spectrum analysis. The color scale indicates excess 
gamma-ray events within an integration radius of 0°.23. The cir- 
cle at the lower left corner shows the size of the point-spread 
function (68% conta inment). The radio contours at 1420 MHz 
HPineault et al.ll997l) are overlaid in white, showing the SNR shell. 
The VERITAS significance contours at 3, 4, 5, and 6t T are shown 
in g reen. The cross marks the position of the pulsar lAbdo et al.l 
1200 8). located 5' ± 2' from the centroid of TeV emission. North is 
up and east is to the left. 

Table 1 

Differential flux measurements of CTA 1 with VERITAS. The 
errors are statistical only. 



Energy Range 


Flux 




Significance 


[TeV] 


[cm-2 s-i TeV" 






0.56 - 1.00 


(1.9 ±0.9) X 10- 


12 


2.3 


1.00 - 1.78 


(7.3 ± 1.7) X 10- 


13 


4.5 


1.78 - 3.16 


(1.2 ±0.4) X 10- 


13 


3.3 


3.16 - 5.62 


(3.4 ± 1.2) X 10- 


14 


3.1 


5.62 - 10.00 


(1.2 ±0.5) X 10- 


14 


2.5 


10.00 - 17.78 


(7.1 ± 2.5) X 10- 


15 


2.8 



radio contours at 1420 MHz, obtained using the Domin- 
ion Radio Astrophysical Observatory (DRAO) Synthesis 
Telescope, and the Effelsberg 100-m telescope, showing 
bright radio arcs visible to the south and east, with an 
incomplete shell in the northwest, possibly due to the 
breakou t of the SNR blast wave in to a medium of lower 
density (jPineault et al.lll993l[l997l ). 

For spectral analysis, the moderate-spectrum cuts 
were used in order to provide the lowest energy thresh- 
old for the analysis. The differential photon spec- 
trum above 500 GeV is shown in Figure [21 with spec- 
tral data points listed in Table [TJ The spectrum is 
generated with the ref lected-region background model 
(jBerge. Funk fc HintonI I2007D with 41 hours 15 min of 
quality-selected data. The spectrum can be fit with a 
power-law of the form dN/dE = No{E/3 TeV)"^, with 
r = 2.2±0.2,tat±0.3sys and Nq = (9.1±1.3stat±1.7sys) x 
10"^"* cm~^ s~^ TeV"^. The integral energy flux above 
1 TeV, = 4.0 X 10~^^ erg cm~^ s~^, corresponds to 
0.2% of the pulsar spin-down luminosity at 1.4 kpc and 
~ 4% of the steady TeV gamma-ray emission from the 
Crab Nebula. 
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E (TeV) 

Figure 2. VERITAS differential gamma-ray spectrum of CTA 1. 
The black butterfly shows the uncertainties of the best-fit 
power-law model. The red line marks the flux predicted by 
IZhane. Jiang fc Lin I ||2009| . Figure 4). The errors arc statistical 
only. 

4.1. Morphology 

Figure [T] shows that the extent of the TeV gamma-ray 
emission region in CTA 1 exceeds the point-spread func- 
tion (PSF; 68% containment radius of the events coming 
from a point source) of VERITAS. In order to estimate 
the extent of the source, an asymmetric two-dimensional 
Gaussian is fit to the acceptance-corrected uncorrelated 
map of excess events binned in 0°.05 bins. Although the 
shape and extent of the emission is likely more complex 
than a simple asymmetric Gaussian, as a first approxima- 
tion, it still provides a statistically reasonable estimate 
of the source extent. Due to the finite resolution of the 
detector, the emission we see is a convolution of the real 
source and the PSF describing the system. Accounting 
for the PSF of the instrument, the resulting la angu- 
lar extent is 0°.30 ± 0°.03 along the semi-major axis and 
0°.24±0°.03 along the semi-minor axis, with an orienta- 
tion angle of 17°. 4 ± 15°. 8 west of north. We note that 
this is a sensitivity-limited measurement. 

The fitted centroid location is 

00'* 06™ 26", +72° 59' 01.0" (J2000), which is 5 
arcmin from PSR J0007+7303. Therefore, the VERI- 
TAS source name is VER J0006+729. The statistical 
uncertainty in the centroid position is 0°.09 in RA and 
0°.04 in declination. The systematic uncertainty in the 
position due to the telescope pointing error is 50" . 

4.2. Archival X-ray Analysis 

Figure [3] shows the exposure-corrected, smoothed 
ROSATFSPC X-ray image of the region around CTA 1. 
The 0.5-2 keV i?OS'^T image reveals a center-filled mor- 
phology, with a faint compact source located at the 
peak of the brightness distribution. The cross in the 
image marks the location of the X-ray point source 
RX J0007.0-1-7303 and the Fermi pulsar J0007-I-7303. 
The pulsar is located close to the center of the ex- 
tended TeV source with 5' offset from the peak of 
the TeV surface brightness. Figure [4] shows the non- 



Right Ascension (J2000) 

Figure 3. ROSAT X-ray image (0.5-2.0 keV) of the SNR CTA 1 
shown in equatorial coordinates. The cross marks the loca- 
tion of the X-ray point source RX J0007.0-f 7303 and the Fermi 
PSR J000 7-I-7303. The SNR s hell is shown by the 1420 MHz radio 
contours HPineault et aI.lll997I V overlaid in white. The VERITAS 
significance contours at 3, 4, 5, and 6(t are shown in green. 



thermal X-ray image from ASCA in the 4-10 keV band 
(|Roberts. Romani fc Kawail [200ll ) . along with the TeV 
contours. While the non-thermal emission and the TeV 
emission are both clearly extended with centroids sepa- 
rated by only 5 arcmin, the limited angular resolution 
and signal-to-noise ratio in the VERITAS data do not 
allow for a more rigorous comparison. 

Figure [5] is the smoothed XMM-Newton image of 
the vicinity of PSR J0007-I-7303 (ObsID 0604940101, 
PI P. Caraveo) showing the X-ray PWN. The inset 
shows the smoothed Chandra image (ObsID 3835, PI 
J. Halpern) revealing a compact nebula and bent jet at- 
tached to the point source, along with diffuse emission 
at larger scales. The Chandra jet is particularly appar - 
ent in the analysis presented by IHalpern et al.l (j2004D , 
where it is estimated that the Chandra point source ac- 
counts for ~ 30% of the flux of RX J0007.0-h7303, with 
the compact nebula plus jet comprising the remaining 
~ 70%. The luminosity of the fainter large-scale emis- 
sion (within r < 4' from the pulsar) is about a factor of 
5-10 larger than that of the compact PWN and pulsar. 
The X-ray spectrum of the point source can be described 
by an absorbed power-law plus blackbody model, with 
a photon index of F = 1.6 ± 0.6. The compact PWN 
spectrum is harder, with a photon index F ~ 1 — 1.3 
(jHalpern et al.l [20M ) . The s pectrum of the larg e-scale 
diffuse emission was fitted bv lCaraveo et al.l (pOlO) with 
a power-law modified by the interstellar absorption. The 
fit gave F — 1.8 ± 0.1. However, the fit quality was fairly 
poor suggesting a more complex spectrum (e.g., a possi- 
ble additional thermal component; see below). 

5. DISCUSSION 

We have discovered spatially extended TeV emission 
from the region of CTA 1. Here we discuss the results in 
the context of the available multiwavelength data. 



Right Ascension (J2000) 

Figure 4. ASCA CIS image (4-10 keV) of the SNR CTA 1, us- 
ing the same field of view as Figures [T] and [3] The position of 
PSR J0007-I-7303 is marked by the cross. The 1420 MHz radio 
contours are shown in white. The VERITAS significance contours 
at 3, 4, 5, and 60" are shown in green. 

5.1. The nature of the TeV source: A PWN scenario 

The good positional match between VER J0006+729 
and PSR J0007+7303 within the CTA 1 SNR makes their 
physical association virtually indisputable, while leaving 
uncertain which component powers the TeV emission. 
The source is unlikely to be related to a gamma-ray bi- 
nary or a background blazar given that the TeV emis- 
sion is extended and non-variable. However, the extent 
of VER J0006-I-729 is also much smaller than that of 
the SNR and, hence, the TeV source does not resem- 
ble several S NRs whose shells ha ve been resolved in TeV 
gamma rays ()Komin et al.|[20lTI ). There still remains a 
possibility that only part of the CTA 1 SNR shell is in- 
teracting with a dense molecular cloud which could cause 
a loc al enhancement in the TeV brightness of the emis- 
sion (jKomin et al.l[20TTI ). However, we do not find any 
evidence for such a cloud at any other wavelengt h, in- 
cluding 60 /im IR or HI (see iPineault et al.l I1993D , and 
the high Galactic latitude of the SNR places it nearly 
250 pc above the Galactic plane, much higher than the 
scale height of molecular clouds. 

VER J0006+729 must then be powered by the young 
PSR J0007-I-7303. Extr apolating the pu lsed emission 
seen by the Fermi-LAT ()Abdo et al.ll2012D . assummg a 
power-law above the break energy of •-^ 4 GeV, gives an 
estimated TeV flux 2-3 orders of magnitude below the 
observed flux. Therefore, the most plausible remaining 
scenario is that the TeV emission originates from the X- 
ray PWN surrounding the pulsar. The PWN consists 
of a bent jet, a compact core, and a large-scale diffuse 
component, as seen in the Chandra and XMM-Newton 
images (see Figure [SJ. The be nding of t he jet (see the 
inset in Figure [5l and also , Halpern et al.l l2004) could be 
caused by the ram pressure of the oncoming medium due 
to the NW-SE pulsar motion or by the interaction with 
a reverse shock propagating NW within the SNR extent. 
Alternatively, a kink instability might be responsible for 
the bending of the jet. Indeed, the Vela pulsar jet shows 
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Figure 5. XMM-Newton 93 ks EPIC/MOSl+2 image (0.5-10 keV; pixel size 4"; smoothed with the r = 12" Gaussian kernel) of 
PSR J0007+7303 and its vicinity showing the X-ray PWN. The inset is the higher resolution 50 ks Chandra ACIS image (0.5-8 keV; pixel 
size 2", smoothed with the r = 4" Gaussian kernel). The arrow in the inset image shows the likely direction of the pulsar proper motion 
(based on the shape of the compact PWN). 



some kink-like shape changes. However, these wiggles 
tend to occur on smaller spatial scales, while globally 
the Vela jet is always (during the last 10 years) seen to 
bend toward one side, likel y due to the pressu re of the on- 
coming ambient material fPavlov "erall[200l . Similarly, 
in the XMM-Newton and Chandra images of CTA 1, ob- 
tained at different epochs, the jet is seen to consistently 
bend in the same direction. 

If the change in the jet morphology is caused by ram 
pressure, then we ca n estimate the press ure from the jet's 
curvature, following iPavlov et al.l ()2003f ). Assuming that 
the jet pressure is dominated by the contribution from 
a magnetic field -B_4 — _B/(10~'* G), for a jet curva- 
ture radius i?curv — 10" and a jet diameter djct ^ 1" 
the pressure estimate is Pram ~ 1 x 10^^^ B^^ erg cm~^. 
Despite the fairly high magnetic field assumed, the es- 
timated pressure is rather low compared to the ambi- 
ent pressure inferred for other yo ung pulsars with X- 
ray P WNe resolved by Chandra (jKargaltsev fc Pavlovl 
12008(1 . Assuming that the ram pressure is caused by 
the pulsar motion through a medium of density n_i = 
n/(10~^ cm""^), one obtains a very modest pulsar speed 

u ~ 90 km s~^ which corresponds to the 

proper motion of just ^ 0.013 n_Y^ B^idi \ arcs yr~^, 
assuming a distance of diA = d/(1.4 kpc). This means 
that over its lifetime r the pulsar should have moved by 

only 3 nZY^ {t/1.3 x 10"* yTs)B^4d^\ arcminutes. This 
distance is much less than the size of the SNR and the 
extent of the TeV source, and it is even less than the ex- 
tent of large-scale X-ray PWN seen in the XMM-Newton 
and ASCA images. We also note that despite being de- 
pendent on several parameters, the above estimate of the 



distance traveled by the pulsar likely represents an upper 
limit. Thus the estimated velocity from the jet-bending 
is inconsistent with the otherwise plausible hypothesis 
that the NW extension of the TeV source (see Figure [1]) 
might be due to aged relativistic electrons left behind by 
the fast moving pulsar. 

It is also possible that the relic PWN has been pushed 
to one side (i.e. NW of the pulsar) by the reverse shock 
that must have arrived from the SE direction. Indeed, 
such a scenario is supported by the overall asymmetry of 
the SNR shell which appears to expand into much lower 
density medium in its NW part and hence the reverse 
shock is not expected to arrive from the NW direction. 
A recent interaction with the reverse shock could pos- 
sibly also explain the bending of the jet while emission 
ahead of the pulsar could be explained by the turbulent 
mixing between the pulsar wind and SN ejecta behind 
the r everse shock (similar to G327. 1-1.1; iTemim et al.l 
120091 ). The latter may contribute some thermal emission 
and explain the poor qual ity of the power-law fit to the 
extended X-ray emission (jCaraveo et al.ll2010l ). A deep 
mapping of the entire SNR by XMM-Newton can test this 
hypothesis by providing a high-S/N spectrum of the faint 
large-scale X-ray emission which should then contain a 
thermal emission com ponent coming from the e jecta (c.f. 
e.g.. Vela X spectrum: iLaMassa. Slane fc de Ja gcr 2008). 

Pulsar wind particles may be transported either by dif- 
fusion or by advection. One could in principle compare 
the two terms, if the bulk flow speed (as a function of 
distance from the pulsar) and the magnetic field struc- 
ture were k nown. The MHD m odels for isotropic pulsar 
winds fe.g.. iKennel fc Coronitil [l984f ) are unlikely to be 
valid on large scales and when mixing due to interaction 
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with the reverse shock is present. However, we can make 
some estimates of the average magnetic field by assuming 
which process is dominant in transporting the particles. 

Assuming that X-ray- and TeV-emitting particles move 
away from the pulsar with similar velocities (i.e. that the 
X-ray- and TeV-emitting regions are co-spatial and that 
the effects of energy-dependent diffusion are negligible), 
and that the synchrotron cooling-time is the dominant 
time-scale, one can crudely estimate the rt iagnetic field 
strength (see, e.g.. lAharonian et al.ll2005b( ). For X-ray 
and TeV gamma-ray emission regions of sizes i?x and R-^ , 
respectively, with Ey^ and being the corresponding 
mean energies of the photons in keV and TeV units, the 
magnetic field is -Bpwn ^ ^(iO{Rx/ Rj)'^ {Ex/ Ej) iiG. For 
the observed Rx/Rj « 0.5, Ex = 5 keV, and E^ = 
5 TeV, the corresponding average magnetic field is ^ 
40 iiG. This is much higher than what is suggested by 
modeling (see below) and also much higher than what is 
seen in other such evolved systems. 

Note that the ratio Rx/R-y ~ 0.5 is likely an under- 
estimate and the ASCA data suggest that it can be a 
factor of 2-3 larger (see Fig. |4|). Indeed, in the ASCA 
images some diffuse ei nission appears to be seen up to 40' 
away from the pulsar (jSlane et al.ll2004 ) . The true extent 
and the nonthermal nature of the faint X-ray emission 
can only be measured in deep observations with XMM- 
Newton. Similarly, the TeV size we quote is a lower limit 
since more the extended portions away from the pulsar 
are likely to be the fainter than our detection threshold. 
The estimate of the average magnetic field should thus 
be taken as a crude order-of-magnitude estimate. 

Another estimate of the magnetic field can be made 
by assuming that diffusion is the dominant transport 
mechanism throughout the nebula. In the simplest case 
of cross field diffusion (the Bohm limit), the diffusion 
constant is given hy D = jmc^ /3eB, where 7 is the 
electron Lorentz factor. Using the relation E^ ^ 7^6, 
where E^ is the mean up-scattered IC photon energy 
and e is the seed photon energy, the diffusion con- 
stant for electrons scattering on the CMB in a mag- 
netic field i?_5 = -Bpwn/(10~^ G) can then be ex- 
pressed as -D = 8.5 X lO'^^Ej^'^Bzl cm^ s^^. Assum- 
ing that particles travel during their characteristic cool- 
ing time of pa 100(1 -|- l-iA B'j^)-^ E'^'"^ kyrs (see 
Ide Jager fc Diannati-Ataj|2009[ ). the diffusion length is 
~ (6Dt^)1/2 _ 13515/^(1 + 14.4B25)-i/2 pc which 

translates into ~ 2Q'd^\ for i3pwn = 5 /^G. (Note that, 
for a given distance to CTA 1, this estimate depends 
only on the magnetic field strength.) This size roughly 
corresponds to the observed extent of the TeV source. 
Although small, such low magnetic field (~ 5 jiG) was 
inferred through the i nultiwavelength spectral mqdeliii g 
for the Vela X plerion ()de Jager. Slane fc LaMassal2008[ ). 
The low Bpwn resulting from the Bohm diffusion esti- 
mates has b een previously noticed for s evera l other relic 
PWN e (e.g. Ide Jager fc Diannati-Atai' (2009), as well as 
bv Uchivama et al. (2009) and Anada et al. ( 201Q) based 
on their analysis of the synchrotron spectra measured by 
Suzaku across the extent of TeV sources HESS J1825-137 
and HESS J1809-193.) 
Furthermore, we can use a dynamical model for 
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Figure 6. Upper : Time evolution of the PWN magnetic field, 
using tlie model of 'Gelfand, Slane & Zliang I ((20091), with parame- 
ters given in Tablc[2l Lower: Time evolution of the modeled SNR 
(blue) and PWN (red) radii. Horizontal dashed lines indicate the 
current values for CTA 1. The vertical green line indicates the age 
at which the measured SNR radius is reached. (See text for model 
description.) 

the evolution of a PWN insid e a non-radiative SNR 
(jGelfand. Slane fc Zhang 1 120091 ) to estimate the physi- 
cal properties of the PWN. The results are illustrated in 
Figure [HI We find that to correctly reproduce the radius 
of the SNR (lower panel, shown in blue) while simulta- 
neously matching the estimated PWN radius (shown in 
red), the current spin-down power, and the total TeV 
fiux, we require an ambient density ng ~ 0.07 cm^'^. 
This is somewhat larger than that estimated from ASCA 
measurements of the thermal X-ray emission (jSlane et al.l 
[1997, 2004), although those measurements were based on 
observations of a small fraction of the SNR shell. 

Also shown in Figure [HI (upper panel) is the time evo- 
lution of the PWN magnetic field for this model. At the 
current age of ^ 10 kyr implied by the SNR radius, the 
PWN magnetic field strength is ~ 6/iG. The recent de- 
crease in the PWN radius, and increase in the magnetic 
field, result from the beginning of the SNR reverse shock 
interaction with the nebula, as suggested by other ar- 
guments presented above. Figure [7] shows the archival 
broadband data for CTA 1 along with the emission pre- 
dicted for the model used in Figure [5] assuming a broken 
power law injection of particles from the pulsar, for which 
a braking index of 3 is assumed. The model parameters 
are summarized in Table [5J 

Radio observations of CTA 1 do not provide 
conclusive evidence for emission from the PWN 
fPincault et al. 19 971). We no te that, in their model- 
ing, Zhang, Jiang fc Lin I ()2009[ ) assumed that the entire 
emission from the SNR was associated with the PWN. 
In fact, the PWN is much fainter. He r e we h ave used the 
1.4 GHz image from iPineault et al.l (|1997[ ) to estimate 
the flux within a 20 arcminute radius around the pulsar, 
and have used this flux as an upper limit for the PWN 
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Figure 7. Broadband emission from CTA 1, along with the syn- 
chrotron and in verse Compton emission from th e PWN predicted 
by the model of IGelfand. alane &: Zhang I II2009I ') with parameters 
given in Table [2l See the text for the derivation of the 1.4 GHz 
radio upper limits shown in blue. The red lines mark the errors on 
the X-ray spectrum measured by ASCA (Slane ct al. 2004). The 
green, das hed lines mark th e errors on the unpulsed Ferrrai-LAT 
spectrum IIAbdo et al.|[20ig) . 

emission. In Figure [71 we have extrapolated this up- 
per hmit to lower frequencies assuming a spectral index 
a — 0.3 (where oc is the flux at the frequency 

v) and to higher frequencies assuming a = 0. These in- 
dex values represent the typical range observed in radio 
PWN, and the associated flux values correspond to con- 
servative upper limits. 

The TeV spectrum predicted bv iZhang. Jiang fc Lin I 
(|2009[) overpredicts the TeV flux and requires an ini- 
tial spin period in the range 30 — 40 ms. In com- 
parison, the model shown in Figure [7] agrees with the 
VERITAS data, within uncertainties, and we conclude 
Pq ^ 1 55 ms. This difference co uld arise from the 
fact that IZhang. Jiang fc Lin I ([20091) assign the full ra- 
dio flux from CTA 1 to the PWN, thus requiring sig- 
nificantly more energy in the electron spectrum. An 
additional difference in the models is that we have cal- 
culated the evolving magnetic field strength based on 
the fraction of spin- down energy inject ed as magnetic 
fiux whereas IZhang. Jiang fc Lin I (|2009f ) assume a time- 
dependent field value which is independent of any other 
system parameters. Our results suggest a break energy 
of ~ 50 GeV with ^ 80% of the spin-down power appear- 
ing in the form of particle flux. Like most such systems, 
CTA 1 is thus a particle-dominated PWN. 

The dashed green curves in Figure [7] represent the 
best fit f or the unpu l sed F ermi-LAT spectrum pub- 
lished in lAbdo etaTI (|2012| ). In the model calcu- 
lation shown, the TeV emission is produced by in- 
verse Compton scattering of photons from the CMB, 
integrated starlight, and infrared emission from local 
dust, following the approxima t e pre scription given by 
iStrong. Moskalenko. fc Reimeii ()2000[) . The model pro- 
duces reasonable agreement with the radio. X-ray, and 
TeV data with a solution that gives approximately the 
correct pulsar spin-down power and characteristic age 
at the current epoch. However, the model is a poor fit 
to the reported Fermi-hAT spectrum. We have consid- 
ered additional photon fields to produce enhanced inverse 
Compton emission at GeV energies, but have been un- 
able to reproduce the published spectral index. We note 



Table 2 

Model parameters for broa dband emission from CTA 1. See 
IGelfand. Slane He Zhangi f|2009i ) for parameter definitions. 



Parameter 




Value 




Input: 




Explosion Energy, -Bgjv 




10^1 erg (fixed) 


Ejecta Mass, M^j 




6.1 M0 


Ambient density, no 




0.068 cm-3 


Initial spin-down, Eq 




7.5 X 10^® erg s~^ 


Spin-down timescale, tq 




3.2 X 10^ yr 


Braking index, n 




3 (fixed) 


VB 




0.2 


ai 




0.5 


012 




2.8 


Break energy, i?f, 




50 GeV 


Age 


Output: 


1.0 X 10* yr 


BpwN 

E 




6.3 ^^G 
4.4 X 10^^ erg 


Tc 




1.3 X lO"* yr 


Po 




155 ms 



that the reported unpulsed GeV emission is quite faint 
and it is in the presence of bright pulsed emission from 
PSR J0007-f 7303. It will be of considerable interest to 
carry out further investigations of this unpulsed emission 
as more Fermi-hAT data are accumulated. 

5.2. Comparison with other relic PWNe 

Figures [8] and [9] show the comparisons of the proper- 
ties of CTA 1 with other PWNe and PWNe candidates 
detected at TeV energies. At the distance of 1.4 kpc, the 
> 1 TeV luminosity of the PWN in CTA 1 is 9.4 x lO^^ 
erg s^^. Figure [8] shows the relative luminosities of 
PWNe in the TeV and X-ray bands, as functions of spin- 
down power and characteristic age (Kargaltsev & Pavlo3 
I2010f l. CTA 1 fits with the picture that TeV PWNe 
are generally found around pulsars with ages ^100 kyrs 
and E > 10^^ erg s~^, although the TeV luminosities 
do not depend on the pulsar age nearly as much as 
the X-ray PWNe luminosities do. Figure |9] shows the 
distance-independent ratio of the TeV to X-ray lumi- 
nosity as a function of characteristic age for a set of 
PWNe or PWN candidates, with CTA 1 marked by the 
red triangle. Additionally, the estimated size of the TeV 
emission region in CTA 1 is consistent with the sizes of 
TeV nebulae around pulsars with ages similar to that of 
PSR J0007-I-7303, although the large errors on the esti- 
mated distances prohibit a definite correlation. A com- 
parison of CTA 1 with the TeV/X-ray PWN population 
supports the PWN origin of the TeV emission. 

6. SUMMARY AND CONCLUSIONS 

VERITAS has detected TeV gamma-ray emission co- 
incident with SNR CTA 1. The emission is ex- 
tended, with a centroid near the Fermi gamma-ray pul- 
sar PSR J00074-7303 and its X-ray pulsar wind neb- 
ula (PWN). The photon spectrum is well described by 
a power-law with differential spectral index of F = 2.2 ± 
0.2stat i 0.3sys and an integral flux above 1 TeV corre- 
sponding to ~ 4% of the steady TeV gamma-ray emission 
from the Crab Nebula. It is unlikely that the TeV emis- 
sion is due to interaction of the CTA 1 shell with a dense 
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Figure 8. Pulsa r spin-down luminosity vs age, from 
[Kargaltscv fc Pavlov! (2010), with CTA 1 point indicated. 
Filled circles: X-ray (red) and TeV (blue) luminosities of PWNe 
or PWN candidates. Larger circle sizes correspond to higher 
luminosities in the correspon ding waveband. Small black dots 
denote ATNF catalog pulsars l IManchester et al.ll2005l )''. 



Figure 9. Plot of the ratio of TeV to X-ray luminosity vs pul- 
sar spin-down age, from Kargaltscv & Pavlov (2010), with CTA 1 
shown by the red triangle. Thick and thin error bars correspond 
to firm and tentative (or questionable) PWN associations, with 
selected objects labeled (see iKargaltsev fc Pavlovi for further de- 
tails). 



molecular cloud, given the lack of evidence for such a 
cloud at other wavelengths (60 /im IR or HI maps). We 
have analyzed archival X-ray data from ROSAT (0.5- 
2 keV) and ASCA (4-10 keV) of the large scale neb- 
ula and XMM-Newton (0.5-10 keV) and Chandra (0.5- 
8 keV) of the region close to the pulsar and find that the 
large scale emission seems to match the TeV morphol- 
ogy. The positional coincidence with the pulsar, small 
extent of the TeV emission, and X-ray/TeV luminosities 
strongly argue for a PWN origin. We have estimated 
the magnetic field strength assuming particle transport 
by either diffusion or by advection. A more detailed dy- 
namical model of the SNR-PWN system suggests a 6 fj,G 
field along with a recent interaction between the PWN 
and the SNR reverse shock. 
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